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Abstract. We present a formal derivation of the key equations governing 
PsJ [ gravitational lensing in arbitrary space-times, starting from the basic properties of 

Jacobi fields and their expressions in terms of the exponential map. A careful analysis 
O . of Jacobi fields and Jacobi classes near the origin of a light beam determines the nature 

of the singular behavior of the optical deformation matrix. We also show that potential 
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original argument presented by Seitz, Schneider & Ehlers (1994). 
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1. Introduction 

X 

Gravitational lensing can be described, starting from very basic principles, in terms of 
the deviations of null geodesies with respect to a "fiducial" ray in a light beam. The 
formal groundworks upon which this description is based were first spelled out by Seitz, 
Schneider and Ehlers (1994) pQ, and their derivation has been widely used ever since 

[2J El H 01]. 

The fundamental objects in this description are the separation vectors which 
determine how a beam of geodesies starting (or ending) at a given point deviates from the 
fiducial. The relevant components of these vectors naturally belong to a 2-dimensional 
space-like surface (the screen) which is orthogonal to the direction of propagation of the 
null fiducial geodesic. And since, by construction, the beam is focused on the reference 
point, the separation vectors are such that £(0) = 0. This description is time-symmetric, 
in the sense that the reference event can be regarded either as the original source of the 
beam (in which case the affine parameter is future-oriented), or as an observation event 
(in which case the affine parameter is past-oriented). 
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The separation vectors are in fact the projection on the screen of Jacobi fields along 
the fiducial ray. A fundamental result in General Relativity is the fact that, to linear 
order in small perturbations around the fiducial geodesic, the Jacobi equation is linear. 
When projected on the screen, that equation leads to £"(A) = T(A) £(A) , where T is 
called the optical tidal matrix, and primes denote derivatives with respect to the affine 
parameter A along the fiducial geodesic. The separation of a null geodesic from the 
fiducial ray at any given value of the affine parameter, £(A), would then be given by 
the action of a linear map (the Jacobi map, [lj) on the velocity of separation of that 
geodesic at the reference point, £'(0) =: 6. Ultimately, these two facts together allow 
us to frame gravitational lensing entirely in terms of the deviations on the screen at 
the reference point. 

We point out that previous demonstrations of these fundamental results have relied 
on a flawed argument which, if taken at face value, would imply that the projections of 
the Jacobi fields on the screen would vanish identically. E.g., the argument presented 
in pQ is the following: since the Jacobi equation is linear, the projection on the screen 
of the Jacobi fields, £(A), can be related to the initial deviation £'(0) = by a linear 
transformation, £(A) = V{\)6. Substituting this expression back in the Jacobi equation, 
they are then able to derive the second-order linear differential equation which governs 
the evolution of the linear operator T>(\). 

They also claim that the linear operator T> obeys the first-order differential equation 
T>' = ST>, where S is called optical deformation matrix and is given in terms of 
optical observables for beams of null geodesies. This would imply, however, because 
the deviation can be written as £' = V 6, that £' = S $, for all values of the affine 
parameter. The problem with this last identity is, of course, that a geodesic with 
£(0) = would necessarily also have to satisfy £'(0) = 0, which would then imply that 
£ = for all values of the affine parameter. The only possible fix for this construction 
is for the optical deformation matrix to be singular at the reference point, in precisely 
the right way to introduce some geodesic deviation £' at that point - but this then 
constitutes an incomplete argument. 

Here we will close the loophole in this argument. First, we will clarify the statement 
that the linearity of the Jacobi equation leads to the Jacobi map. Although our argument 
also relies on the linearity of the Jacobi equation, we present an explicit construction 
of the relation £(A) = T>{\)6 where T>(\) is directly determined from the exponential 
map. Second, by employing the notion of Jacobi tensors we will show how to obtain 
a second-order linear differential equation for T> which depends only on the Ricci and 
Weyl curvature tensors. 

In our demonstration we will use a few basic results from Lorentzian geometry, 
in particular the fact that basis vectors along null geodesies cannot be orthonormal, 
and decomposition of vectors and tensor in terms of these basis are not defined in the 
classical sense. We follow [7] in our construction of such a basis, over which separation 
vectors can be expressed. For this basis, we need to introduce the notions of the quotient 
by an equivalence relation. This then allows us to employ, in Section HI Jacobi classes 
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and Jacobi tensors [SJ, which turn out to be key to our argument. 

In Section [5] we construct the Jacobi map, and obtain the differential equation that 
is satisfied by the operator T>{\). After writing the basic equations that govern the 
lensing problem in terms of the exponential map, we reexamine the argument presented 
by pp. Finally, we present some comments about conjugate points and critical behavior 
of beams of null geodesies. 

Our conventions are as follows. Space-time is assumed to be a Lorenzian manifold 
(M, g), and the tangent space to a point p G M will be denoted by T P M. Our metric has 
signature — , +, +, +, and contractions with it are denoted as (.,.), i.e., (x, y) = g a px a y^ '. 
Greek indices range from to 3 and latin indices from 1 to 3, except the indices i and 
j, which can only take the values 1 or 2. The 2x2 identity matrix will be denoted by 
1. The Riemann and Ricci tensors are denoted by R( . , . ) . and Ric( . , . ), respectively. 



2. The exponential map and Jacobi fields 

Let 7 : [0, 1] -* M be a null geodesic, so, by definition, § (g) = and (g, g) = 0. 
For given initial conditions 7(0) = p, and 7'(0) = k G T p M, the solution of the geodesic 
equation leads to a curve 7(A,p, k). If the solution 7 is defined, the affine parameter 
and the initial velocity can be scaled in such a way that 7(1, p, k) G M. The point 
7(1, p, k) G M is called the image of k under the exponential map, for p fixed. In our 
notation, the exponential map exp p : T P M — > M is defined for a given k G T P M as the 
point obtained along the unique geodesic in M with 7(0, p, k) = p and 7'(0,p, k) = k, 
after displacement by a length of 1 along this curve, i. e., exp p {k) := 7(1, p, k) [8J. An 
illustration of the action of the exponential map is shown in figure [TJ 

Determining the explicit form of the exponential map is, in general, a very difficult 
task, and relies on the knowledge of all geodesies for a given space-time. Even without 
explicit formulas for the exponential map, however, much can be learned from its 
general properties. The geodesic homogeneity lemma, for example, establishes the 
aforementioned scaling relation between the affine parameter and the velocity, i. e., 
7(1, p, k) = j(l/X,p, \k), and allows us to calculate the derivative of the exponential 
map at A = 0: for any k G T p M, 

-^[exp p (A£;)] A= o = ^[7(1, p, Xk)} x = = ^t[7(A,P, k)] x=0 = k . (1) 

The link between exponential map and variations of curves in the manifold is largely 
employed in the context of variational calculus [9], and can be established as follows: 
let k(s), s G (—£,£) C 1, be a curve in T p M such that A;(0) = k, dk/ds(0) := w. 
Then f(X,s) := 7(1,2?, Xk(s)) = exp p [AA;(s)] defines a parametrized surface in M. For 
s = we obtain the curve 7(A,p, k), and other values of the parameter s generate 
variations of this original curve. A variational vector field J along 7 constructed as 
J(A) = |£(A,s = 0) is such that J(0) = and J'(0) := J^(A,s = 0) = w. Such a 
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Figure 1. Action of the exponential map. Given a vector in T p M, its image under the 
exponential map is the point in M obtained after a displacement of length 1 along the 
(unique) geodesic starting at p with a that vector as velocity. We show the image of the 
point p under the exponential map for three different initial velocities u,k,v G T p M. 



vector field is known as Jacobi field and satisfies the Jacobi equation [7]: 

^+ J R( 7 '(A),J(A)) 7 , (A) = 0. (2) 

Despite being broadly interpreted as geodesic deviation equation, the Jacobi 
equation have other solutions, depending on the initial conditions. For instance: if 7(A) 
is a geodesic, then 7'(A) is a Jacobi field along 7, because 7'(A) is parallel-transported 
along itself, and because the Riemann curvature is anti-symmetric with respect to its 
two first arguments. Another example of a non-trivial Jacobi field along the curve 7(A) 
is A 7 '(A). 

Describing gravitational lensing, however, requires to take into account 
deformations of beams of null geodesies that begin at a given specific point in space-time 
and, therefore, conduce us to consider Jacobi fields that obey initial conditions J(0) = 
and J'(0) = w, where w is the initial velocity of separation between a given geodesic 
in the beam and the fiducial one. A fundamental property of these Jacobi fields, which 
lies at the heart of our discussion, is that the unique Jacobi field along 7 with J(0) = 
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and J'(0) = w is given by - see, e.g., Prop. 6 in chapter 8 of [10], or Prop 10.16 of [SJ: 

J(A) = A(cfexp p ) Afc w. (3) 

To make more clear the notation employed in [3] and show that it defines, in fact, a 

Jacobi field, we can explicitly calculate: 

df d 
J{a) = -r^(\s)\s=o = ^exp p (AA;(s))| s=0 

( dk \ 

= (dexp p ) Xk (o) ( = °) J = K dex V P )\kW , (4) 

where we have used the chain rule to establish the third equality. Expressing 
the exponential map in terms of its definition, we can write X(dexp p )x k w = 
Jj7(l, P, Xk(s))\ s=0 . The differential of the exponential map (dexp p )\ k is the linear map 

(dexp p )xk : T Afc (T p M) -)• T cxPp ( Afc )M 

whose matrix form is given by the directional derivatives of exp p in the direction of the 
basis vectors of T\k(T p M), evaluated at Xk. In general, for v G T P M, the tangent space 
T v (T p M) is defined as: 

T v (T p M) := {§ w : R ->• T p M} , = v + sw v, w G T p M . 

Therefore, the space T V (T P M) is isomorphic to T P M, so for our purposes the two can 
be identified. Consequently, (cfexp p ) Afc will be understood as a linear map between the 
spaces T P M and T exPp (Afc)^f- In figure [2] we illustrate the action of the differential of the 
exponential map on vectors in the tangent space of a point p G 7. 

To determine the behavior of the differential of the exponential map near the origin, 
we shall rewrite HI inserting an intermediate equality (application of the chain rule) 
between the first and the last one: 

-^[exp p (\k)}x=o = (cfexpp)o k = k. 

In other words, the differential of the exponential map calculated at the origin equals to 
the identity matrix. This can be rephrased as the statement that the exponential map 
defines a diffeomorphism of a neighborhood of the origin of T p M into an open subset of 
M. This property also allows us to verify that J given in [3] does satisfy J'(0) = w: 

J'{0) = ^[X(dexp p )xk w] x =o = (dexp p )x k w\x= + A-^[(cfexp p ) A/c w] x =o = w. (5) 
3. Construction of screen vectors 

We now turn our attention to the Jacobi fields that can be expressed as in [31 and the 
way they can be decomposed in terms of a vectors basis along the null geodesic 7. Let's 
consider the vector basis along a null geodesic as introduced in Sec. 4.2 of [7]: since 
p = 7(0), we define a non-orthonormal basis {Eq, E±, E2, E3} at T P M by requiring that 
i? = 7'(0). By construction, (E ,E ) = 0, and we take E 3 to be another null vector 
that obeys {E Q , E 3 ) = — 1. The remaining vectors E\, E 2 are then space-like and can be 
normalized. To extend this basis along 7 we simply parallel-transport the basis vectors. 
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Figure 2. The action of (dexp p )k- For a given w € T p M, the action of (cfexp p )fc gives a 
vector in the tangent space of the point exp p (fc) € M. The image of w under (dexp p )k 
is interpreted as the separation, at A = 1, of a fiducial null geodesic from another null 
geodesic starting at the same point but with an initial velocity of separation w. Here 
we identify T k (T p M) with T p M. 

With the help of this basis along 7, we can describe the separation of a bundle of 
geodesies that emerge from p with any velocities. We will be primarily concerned with 
the deformation of light beams with respect to the fiducial ray 7, and will describe these 
other geodesies in terms of the basis defined in the tangent space of all points along the 
fiducial curve. Since the separation of geodesies is described by the Jacobi fields, and 
all Jacobi fields that satisfy J(0) = 0, J'(0) = w have the form given in El then all the 
dependence on the separation between the geodesies is codified in J'(0). If we decompose 
J'(0) in the basis {E , E 1 , E 2 , E 3 }, the component J£(0) is such that (4(0), V(0)) ^ 0. 
Hence, it follows from Gauss' Lemma that (^(A), 7'(A)) 7^ for all A, where J3 is the 
component of J in the direction E 3 . This component brings no information whatsoever 
about the problem we are interested in, so we shall only consider the vector space that 
results from the inverse image of by the application (.,7' (A)}. With a notational 
abuse, we shall call this space T P M, even if p was previously assumed to be fixed. 

There is, however, still one null Jacobi field remaining in T p M: as we have already 
remarked, 7'(A) is itself a Jacobi field along 7, but it does not describe any spatial 
separation between geodesies. In order to eliminate this component, consider two 
vectors W and V in T p M, and define W ~ V if W — V £ span(7'). It is easy to 
show that this is an equivalence relation in T P M. The set of equivalence classes in T P M 
by the relation ~ inherit the structure of vector space from T p M. We shall denote 
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Figure 3. Intuitive idea behind the construction of the screen. Any vector to which a 
scalar multiple of k is added, will represent the same spatial separation from k. They 
should therefore be not distinguished for the sake of our problem and, for that reason, 
they are all identified inside a class of equivalence. The set of linearly independent 
equivalence classes forms T p M/span(7'). 

T p M/span(7') the set of equivalence classes of T P M by the relation ~. Given some 
V G T p M, [V] G T p M/span(7') will denote the class corresponding to it. The same 
vector [V] corresponds to the family of vectors V + aj 1 G T P M, a G R, as illustrated in 
figure [3j The two-dimensional space T p M/span(7') we have just constructed is called 
screen, and its elements are called screen vectors. 

The result of the construction above is a two-dimensional vector space associated 
to every point along 7, in terms of which separations of geodesies in the beam can be 
described. All vectors in this space are space-like, and correspond to the naive idea 
of projection to a space "orthogonal" to 7'(A) - after taking into account that 7'(A) is 
orthogonal to itself. 

4. Jacobi classes, Jacobi tensors and optical scalars 

We have just constructed a vectorial space to which the objects we want to describe will 
be restricted. We should ask how the space-time geometry is expressed when restricted 
to this space. Firstly, we note that, given U, V G T P M, then, since T P M has the structure 
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of f p M/span(7') © span^'), U = [U] + ak and V = [V] + (3k for same a, (3 £ R. 
Therefore (U,V) = ([U] + ak, [V] + /3k) = ([U], [V]}. Similarly, Vy[E7] = [VyE/]. Also 

[^(c/^OY] = [R{[uii)i + ocR{i,i)i] = [R{[u],i)i] ■= [R\([u\>i)i, due to the 

symmetries of Riemann tensor. 

We shall now consider the Jacobi fields restricted to the quotient space or, Jacobi 
classes, as these restrictions are usually called. In general, Jacobi classes are smooth 
vector fields in the quotient space satisfying the relation: 

j"+[R}(j,7'W=W]- (6) 
It can be shown that, given a Jacobi class j £ T p M/span(7'), there exists a Jacobi field 
J in T P M such that [J] = j, and, conversely, if J is a Jacobi field in T P M, j = [J] is a 
Jacobi class in T p M/span(7') [8]. 

We now introduce objects called Jacobi tensors [8], which are of fundamental 
importance in the discussion about optical properties of light beams, and in terms of 
which the optical scalars will be defined. Let A be a bilinear form A : T p M/span(7') — > 
T p M/span(y). For any W £ f p M, [R]A([W]) := [R]{A[W],i)i . We shall say that A 
is a Jacobi tensor if: 

A" + [R]A = , (7) 
subject to convenient initial conditions, and 

ker{A(A)} n ker{A'(A)} = [ 7 '(A)] . (8) 

If the bilinear form A satisfies d then vector fields constructed by its action 
on vector basis E t will be Jacobi classes. The condition [8] eliminates trivial Jacobi 
classes from consideration: if ji = AEi are Jacobi classes along 7, since Ei are parallel 
propagated, j[ = A'Ei give their (covariant) derivatives. If at some point along 7 the 
same linear combination of E^s belong to the kernel of both A and A', then the action 
of A on this linear combination will give rise to a trivial Jacobi class. To see that this 
is the case, we should recall that the dimension of the kernel of A is only non-null at 
conjugate points [8], and conclude, by taking a point A* conjugate to A = along 7 for 
which j'(X*) = [7'], that the only possible solution for [6] in this case does not represent 
separation of curves in a beam. 

The condition |8] has an important consequence: objects defined by A' A -1 may be 
singular either because A is singular, or because A' is singular, but never because both 
are singular at the same point. 

The optical scalars expansion, vorticity and shear of a geodesic beam are defined in 
terms of the Jacobi tensors, namely, through combinations of the object B =: A' A^ 1 . 
Explicitly, we have [8]: 

(i) expansion 

= Tr{B) , (9) 

(ii) vorticity 

w = 1m(B), (10) 



The Jacobi map for gravitational lensing: the role of the exponential map 



9 



(iii) shear 

a = Re(£>)-|l. (ii) 
In terms of these objects the Raychaudhuri equation reads: 

0' = -Ric( 7 ', i) - Tr(u, 2 ) - Tr(<x 2 ) - ^ . (12) 

Jacobi tensors are, then, key objects to understand and describe deformations of 
light beams. Next we will address the explicit construction of a Jacobi tensor that 
satisfies the physical restriction: since we are interested in gravitational lensing, all 
geodesies in the beam start at the same point and have different separation velocities 
with respect to a fiducial ray. 



5. The Jacobi map 

In what follows we will explicitly built the Jacobi map. The argument of Ref. [1] is that 
the linearity of the Jacobi equation implies that the Jacobi map should hold - although 
this is in fact true, we will see that it is far from trivial. 
By definition, in T p M/span(7'), 

ji(X) = ([J(A)], [Ei(X)]) = (J{\),Ei(\)) = ((\dexp p Uw,Ei(\)) , 

where w = J'(0) and i = 1,2. This initial velocity can be expressed as w = 
9iEi(0) + 9 2 E 2 (0) + aE and, since ((Xdexp p )\kE , Ei(X)) = 0, there is no loss of 
generality in taking w to be restricted to T p M/span(7'). This just means that we 
only have to determine the space-like components of the initial velocity of dispersion 
of geodesies in order to determine their future evolution. The tangential component of 
the geodesies may also deserve attention, for instance in cosmological contexts where 
redshifts or the Sachs- Wolfe [TJJ effect can take place, but in those cases the corrections 
can be treated separately. 

Taking the initial spread velocity restricted to T p M/span(7'), then ([A (dexp p )\k ■])[■]) 
constitutes a bilinear form in T p M/span(7'), and therefore it admits a matrix represen- 
tation: 

h\ = ( (X(dsxp p ) }Je E 1 (0),E 1 (X)) (A (dexjp p )\kEi(0), E 2 (X)) \ f B x 
j 2 J \(X{deicp p ) xk E 2 {0),E 1 {X)} (X(dex Pp ) Xk E 2 (0),E 2 (X)) ) { 9 2 

IT31 expresses the separation of null geodesies in terms of the initial velocity of 
dispersion of the beam. If we denote the 2x2 matrix of [13] by T>(X), even if it depends 
on the point p (A = 0) and on the fiducial geodesic 7, we obtain [j: 

j=X?(A)0. (14) 

This is, in fact, what is usually called Jacobi map in the gravitational lensing 
literature []]. Using the explicit form given in [13], it is easy to see, after [3] [1] and EJ that 
the matrix T> satisfies the initial conditions T>(0) = and T>'(0) = 1. 

I Ref. PQ calls £ what we are calling j . 
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6. A differential equation for V 

We will now show that the matrix T> which appears in the Jacobi map, [HI satisfies the 
differential equation [7] This can be verified with the help of the Jacobi equation [2j 

Let's define e = [E±\ + i[E 2 ], J = ji + ij 2 , and let e and J be their complex 
conjugates. Since for a given Jacobi field there is a unique Jacobi class associated with 
it, 

(e, [J]"} = - (e, [R(J,i)i]) = [RW,i)i) 

- j±(e, [R](en'W) = -nj - TJ (15) 

where 1Z and T are given in terms of the Riemann curvature. Writing k = 7', it can be 
shown [TJ] that these objects are given by 

U^{e,[R]{e,k)k) = X -R^k v , (16) 

and 

T = i(e, [R](e, k)k) = \e a e?C a ^ v k»k v , (17) 

where R pu and C afl/ 3 U are the components of the Ricci and Weyl tensor, respectively. 
Substituting [TJ] in [15] we obtain, indeed, [7J Explicitly, 

V"{\) + T{\)V{\) = (18) 

where T is the matrix form of the endomorphism [R]( . acting on T p M/span(7'), 

and is usually called optical tidal matrix. In terms of TZ and J 7 , T reads: 

r- ( n + Re:F lmjr ) (i9) 

As already remarked, [18] must be subjected to the initial conditions ^(0) = and 
V'(0) = 1. 

7. The condition on the kernels 

In order to verify that T> is in fact a Jacobi tensor, we must guarantee that the condition 
expressed by [8] is satisfied by the bilinear form constructed in [13] for all values of A for 
which 7 is defined. Since M holds for A = because of the initial conditions, we shall 
consider the case A 7^ 0. Keeping in mind that the vectors £7, (A) are parallel-transported 
along the fiducial ray, taking the derivative T>'(\) leads to 



V>{\) 




((dexp^xkE^O), £i(A)> ((dexp p ) AjfcJ B 1 (0), E 2 {\)) 
((dex Pp ) Xk E 2 (0), E 1 (\)) ((dex Pp ) Xk E 2 (0), E 2 {\)) 

((dexp^^O), ^(A)) ((dexp^xkE^O), E 2 (X)) 
((dex Pp ) Xk E 2 (0), E 1 (X)) ((dex Pp ) Xk E 2 (0), E 2 {\)) 



.(20) 
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In order to determine if there are common elements in the kernels of T> and V, we 
shall study if elements in the kernel of T> may also be in the kernel of T>' . Given its 
explicit form of [131 an d the fact that Jacobi fields expressed by [3] are linearly independent 
if the initial velocities are linearly independent, ker(D(A)) has non-trivial elements only 
at points which are conjugate points to A = along 7. 

At conjugate points, the first term in the right hand side of [20] is singular. The 
second term in the right hand side of [201 however, must not be singular for the same 
Jacobi class because, if it were, then the Jacobi class thus obtained would be the trivial 
one. We must, therefore, make the derivative term in the right-hand side of [20] to be 
singular because of the second Jacobi class, i.e., the second Jacobi class must reach an 
extremum exactly at the same value of the parameter where the first class vanishes. 
Although this is not sufficient to cause a violation of the condition expressed by [HI if we 
could have a sequence of conjugate points along the two Jacobi classes, as illustrated 
in figure H] (i. e., between any two conjugate points of a class, there is a maximum 
of the other class), and if we were able to make all conjugate points accumulate in a 
neighborhood of some point in the interval, then there could be a way of taking this limit 
such that the condition [8] would fail in some point of that neighborhood. The Morse 
Index Theorem, however, states that the set of conjugate points along a null geodesic 
is a finite set [8]. Consequently, an accumulation of conjugate points cannot take place, 
and hence the limit depicted in figure H] cannot exist. We conclude, then, that E] holds 
for the bilinear form V, and therefore it determines a Jacobi tensor. 

8. The behavior near the origin and the sweep method 

We now return to the argument presented in [1] to show that, despite the loophole 
in their argument leading to an equation equivalent to [THl this is of no consequence. 
Specifically, we will show that a matrix S such that T>' = ST> and j' = Sj can be defined, 
but that it cannot be taken as a starting point for obtaining [TBI 

Suppose that there exists a matrix S such that T>'(\) = S(\)V(\), and that we 
restrict the domain of A in such a way that there are no conjugate points to A = along 
7 . Then, from [TBI it follows that S' + S 2 = T. This equation is known as matrix Riccati 
equation, and corresponds to the usual association of Riccati's equation to second-order 
homogeneous differential equations. 

The initial conditions for this matrix Riccati equation should be carefully 
considered. If both X>(0) = 0, V'(0) = 1 and V'{\) = S(X)V(X) hold simultaneously, 
then 5(0) must not be limited. Consequently, one cannot guarantee the existence of 
solutions for this first-order linear differential equation if A is in the neighborhood of the 
origin (A = 0), and in that case one cannot assure a solution for the associated Riccati 
equation either. 

To avoid problems with the solution around A = 0, we take A > e > and impose 
initial conditions T>(e) and T>'(e), as illustrated in figure [5] The equation V = ST> will 
then admit solutions for A > s, if 5(A) is bounded for A > e. If some further conditions 
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Figure 4. Two independent Jacobi classes along a fiducial ray. Let's say that after 
a given value of the afhne parameter along the fiducial ray, there is a sequence of 
conjugate points such that, between two conjugate points of the same Jacobi class, 
the second class passes through a point of maximum separation. A collapse of such 
a sequence of conjugate points in a neighborhood, depicted by the thick line in the 
diagram in the right, could give rise to a violation of condition [8j This, however, is 
not possible because the set of conjugate points along 7 is finite, as states the Morse 
Index Theorem. 

V(X) and V(X), A > e, are also given, then one can solve for V = SV, subjected to 
these conditions, and check whether the solutions match those obtained by imposing 
initial conditions at A = e. If the two solutions match, the initial conditions given at e 
and A are said consistent, and the solution T> will also be a solution of HE] in the same 
range of parameters. This method of obtaining solutions is known as "sweep method" 

The sweep method would work well if geodesies in the beam did not cross each 
other - in other terms, if they formed a congruence. If the beam emerges from the same 
point, however, initial conditions at that point cannot be used for the "forward sweep" 
step of the method, because one cannot guarantee the existence of solutions for the 
differential equation in this case. In spite of that, if one knows, by any other method, a 
solution for [THJ then the behavior of this solution near A = can be investigated. We 
presented explicitly how D is given in terms of the exponential map in [131 an d therefore 
we can investigate the properties of T>(e), e ~ 0. 

If one recalls that the entries of [13] are projections of Jacobi fields in base vectors, 
and that Jacobi fields in the neighborhood of the origin behave like e + Oe 3 , we see that 
for e close to zero, T>(e) ~ el. Then, from [201 we obtain that T>'(e) ~ 1. Hence, an 
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Figure 5. Beam of null geodesies starting at a point. At A = e the beam behaves like 
a congruence, and an equation like T>'(\) = S(X)T>(\) can be solved for given initial 
conditions, if 5(e) is bounded. If we solve the same equation for initial conditions at 
A = A, we can see if the solutions match and, if they do, then we can say that it is also 
a solution of the second-order equation V" + TV = 0. 



equation like V = ST> would only make sense if near A = the matrix 5 had a behavior 
like 5(e) « (l/e)l. 

We remark that, despite the fact that the equation V = ST> may not be inconsistent 
with the initial conditions T>(0) = and T>'(0) = 1, if 5 diverges in the prescribed way 
near the origin, the initial conditions for the Riccati equation cannot be provided at the 
origin, and the sweep forward cannot be employed. In other words, a general solution 
for [TH] must be known in order to verify that V = SV> is consistent near the origin - 
but this last equation could never be used to derive or to generate solutions to [18] if the 
origin (or conjugate points) are in the domain of the solution. 

9. Conjugate points and critical behavior 

The set of conjugate points along the geodesic 7 is not degenerate when the restriction 
to the quotient space is taken: it is in fact completely preserved, because the Jacobi 
field in the direction of 7' (i.e., 7' itself) does not vanish. All possible non-trivial Jacobi 
fields that satisfy J(0) = 0, and vanish at some other point 7(Ai), correspond to Jacobi 
classes satisfying [J(0)] = [7'(0)] and [J(Ai)] = [7'(Ai)]. The conjugate points along 7 
correspond to the critical points of the matrix T), and induce optical critical behavior 
which is codified by the diverging expansion given by [91 

Because T p M/span(7') is bi- dimensional, the maximal multiplicity of any conjugate 
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point is two. After a conjugate point of multiplicity one, the source's image will appear 
inverted - something that does not occur after conjugate points of multiplicity two (or 
two conjugate points of multiplicity one). 

The presence of conjugate points along 7 is much more significant in aspects other 
than classifying an image's orientation. [T4] showed that conjugate points are necessary 
to form multiple images, and due to the connections between the existence of conjugate 
points along geodesies and the energy conditions [7] , multiple images will always occur in 
a large class of space-times, as long as geodesies can be enough extended. As conjugate 
points are critical points of the energy functional, they also play an important role in 
the formulation of Morse Theory applied to light rays, in particular in the formulation 
of theorems on the possible number of images in lensing configurations [15J. 

10. Discussion 

The description of gravitational lensing presented here is central in many applications, 
especially in cosmology ((TJ El [3J H]), when [TBI can be solved perturbatively to generate 
a lens map. Equations with the same general structure are also fundamental in the 
discussion of singularities in the Penrose limit [H] and can have other applications. 
Writing these equations in a rigorous way has allowed us to determine their explicit 
form in terms of the exponential map, beyond generating a clear interpretation for the 
Jacobi map and validate its centrality in the description of gravitational lensing. 

The limitations of the derivation presented in this work should also be clarified. 
First, the Jacobi equation itself is only approximative in first order to the problem 
of geodesic spread. Corrections may be included, such as non-linearities with respect 
to the velocity dispersion [17]. Second, despite these attempts to enlarge the limits of 
application of generalizations of the Jacobi equation, there are obstacles that may not be 
overpassed in this formulation, such as configurations in which cut points j8] exist, and 
a more careful geometrical analysis is required. One such example is lensing produced 
by a string, as introduced by [18], which induces a flat space-time that nevertheless 
generates multiple imaging. In full generality, gravitational lensing phenomena involve 
conjugate points and even cut points, and we can expect that an observer's past light 
cone will eventually reveal this complex geometry. 
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